Self-energy and lifetime of Shockley and image states on Cu(100) and Cu(lll): 
Beyond the GW approximation of many-body theory 
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We report many-body calculations of the self-energy and lifetime of Shockley and image states 
on the (100) and (111) surfaces of Cu that go beyond the GW approximation of many-body theory. 
The self-energy is computed in the framework of the GWY approximation by including short-range 
exchange-correlation (XC) effects both in the screened interaction W (beyond the random-phase ap- 
proximation) and in the expansion of the self-energy in terms of W (beyond the GW approximation). 
Exchange-correlation effects are described within time-dependent density-functional theory from the 
knowledge of an adiabatic nonlocal XC kernel that goes beyond the local-density approximation. 
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I. INTRODUCTION 

At metal surfaces there exist specific electronic states 
not present in the bulk, which can be classified as intrin- 
sic (crystal-induced) surface states 1 - and image-potential 
(Rydberg-like) states* 2 -' 3 Intrinsic surface states are origi- 
nated by the symmetry breaking at the surface, they have 
their maximum near the surface, and they are classified 
as Tammi and Shockley 5 - states; in particular, intrinsic 
Shockley surface states typically occur in the gap of free- 
electron-like s,p bands near the Fermi level j&Zi Image- 
potential states appear as a result of the self-interaction 
that an electron near the surface suffers from the polar- 
ization charge it induces at the surface, and they occur 
in the vacuum region of metal surfaces with a band gap 
near the vacuum levels 

Figure Q] illustrates Shockley and image-potential 
states in the gap of the TL projected band structure of 
the (100) and (111) surfaces of Cu. If an electron or hole 
is added to the solid at one of these states, inelastic cou- 
pling of the excited quasiparticle with the crystal, which 
can be experimentally observed through a variety of 



spectroscopies,' 
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may occur through electron- 



electron (e-e) and electron-phonon (e-ph) scattering. The 
decay rate due to the e-ph interaction, which is relatively 
important only in the case of excited Shockley holes near 
the Fermi level, has been investigated recently by using 
the Eliashberg function^ Accurate many-body calcula- 
tions of the decay rate due to the e-e interaction were first 
carried out for image states on the (100) and (111) sur- 
faces of Cuj 16 : 17 Since then, many-body calculations of 
the e-e decay have been reported for a variety of simple, 
noble and transition metals i 18 ' 19 : 20 : 21 Nevertheless, ex- 
isting calculations have been typically performed within 
the G°W° approximation of many-body theor y 22 : 23 : 24 : 25 
with no inclusion of exchange and correlation (XC) ef- 
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FIG. 1: The TL projected bulk band structure (shaded areas) 
of the (111) and (100) surfaces of Cu. The solid lines represent 
Shockley (n = 0) and image-potential (n — 1) surface-state 
bands. 



fects. Exceptions are (i) a calculation of the e-e decay 
rate of image states on the (100) and (111) surfaces of Cu 
that incorporates XC effects in an adiabatic local-density 
approximation (ALDA)! 7 - and (ii) an approximate evalu- 
ation of the lifetime of Shockley states in the noble metals 
that incorporates the exchange contribution to the self- 

20 

energy^ 

In this paper, we report extensive calculations of the 
screened interaction, the self-energy, and the e-e inelas- 
tic lifetime of Shockley and image states on the (100) 
and (111) surfaces of Cu that go beyond the G°W° ap- 
proximation. Short-range XC effects are incorporated 
both in the description of the dynamical screening of the 
many-electron system [we go beyond the random-phase 
approximation (RPA) in the evaluation of the screened 
interaction W] and in the expansion of the electron self- 



2 



energy in terms of W [we go beyond the GW approxi- 
mation]. This is the GWT approximation of many-body 
theory , 26 ! 27 which treats on the same footing XC effects 
between pairs of electrons within the Fermi sea (screen- 
ing electrons) and between the excited electron and the 
Fermi sea. 

Mahan and Sernelius^ showed that the inclusion, 
within the GWT approximation, of the same vertex func- 
tion in the screened interaction and the numerator of 
the self-energy yields results for the band-width of a ho- 
mogeneous electron gas very similar to those obtained 
in the G°W° approximation, due to a large cancellation 
of vertex corrections. Large cancellations were also ob- 
served to occur in the decay rates of imaged and bulk 28 
states in the noble metals, by incorporating XC effects in 
the ALDA. In the decay of low-energy bulk states below 
the vacuum level energy transfers Tiu are well below the 
Fermi energy and momentum transfers Tiq are typically 
smaller than 2hqF, Qf being the magnitude of the Fermi 
wave vector, so that one can safely assume that both q 
and u> are small and XC effects can, therefore, be incor- 
porated in the ALDA. However, in the case of Shockley 
and image states the ALDA might lead to spurious re- 
sults, due to the presence of small local values of the 
Fermi wave vector in a region where the electron density 
is small. Hence, here we use an adiabatic nonlocal XC 
kernel that accurately describes XC effects in the limit 
of a homogeneous electron gas of arbitrary density and 
which has been succesful in the description of the XC 
contribution to the jellium surface energy^ 

It has been argued in the past that a realistic first- 
principles description of the electronic band structure is 
of key importance in the determination of the inelastic 
lifetime of bulk electronic states in the noble metals^ 
The main conclusion drawn in Ref. was that in the 
case of the noble metals deviations from electron dynam- 
ics in a free gas of sp electrons mainly originate in the 
participation of d electrons in the screening of electron- 
electron interactions. The role of occupied d bands in 
the dynamics of excited surface-state electrons and holes 
on silver surfaces was later investigated via a polarizable 
medium giving rise to additional screening^ and it was 
concluded that d electrons do not participate significantly 
in the screening of the interaction between surface states 
(which are located near the surface) and the Fermi gas 
of the solidi^ 2 . 

In order to investigate the dynamics of Shockley and 
image states on Cu surfaces, we assume that the density 
of valence electrons in the solid varies only along the z 
axis, which is taken to be normal to the surface. Hence, 
our calculations start by solving the single-particle time- 
independent Schrodingcr equation of electrons moving 
in a physically motivated one-dimensional (ID) model 
potential that is known to correcly reproduce the be- 
haviour of sp valence states and accurately describes, in 
particular, the projected bulk band gap and the bind- 
ing energy of the Shockley and the first image stated 
The eigenfunctions and eigenvalues of such an effective 



single-particle hamiltonian are then used to compute the 
screened interaction, the GWT self-energy, and the e-e 
decay rates of Shockley and image states. For compar- 
ison, we also compute G°W°, G°W, and GW°T decay 
rates, with no inclusion of XC effects, with inclusion of 
XC effects beyond the RPA in the screened interaction W 
alone, and with inclusion of XC effects beyond the G°W° 
in the expansion of the electron self-energy in terms of the 
RPA screened interaction W a , respectively. Our results 
indicate that (i) although the use of the ALDA leads to 
spurious results for the screened interaction, a more re- 
alistic adiabatic nonlocal description of XC effects yields 
inelastic lifetimes of Shockley and image states that esen- 
tially coincide with those obtained in the ALDA, and (ii) 
the overall effect of short-range XC is small and GWT 
linewidths are close to their G°W° counterparts, as oc- 
curs in the case of low-energy bulk states.— 

The paper is organized as follows. Explicit expres- 
sions for the e-e decay rate of surface-state electrons and 
holes at solid surfaces are derived in Sec. II, in the GWT 
approximation of many-body theory. The results of nu- 
merical calculations of the screened interaction, the self- 
energy, and the decay rate of Shockley and image states 
on the (100) and (111) surfaces of Cu are presented in 
Sec. III. The summary and conclusions are given in Sec. 
IV. Unless stated otherwise, atomic units (a.u.) are used 
throughout, i.e., e 2 = fi = m e = 1. 

II. THEORY 

Let us consider an arbitrary many-electron system of 
density no(r). In the framework of many-body theory, 
the decay rate (or reciprocal lifetime) of a quasiparti- 
cle (electron or hole) that has been added in the single- 
particle state <^>i(r) of energy Si is obtained as the projec- 
tion of the imaginary part of the self-energy £(r,r';£i) 
over the quasiparticle-state itself 2 ^ 

V 1 =TlJdrJ <irV>*(r)Im£(r, r'; ^(r')* (1) 

where the =F sign in front of the integral should be taken 
to be minus or plus depending on whether the quasi- 
particle is an electron (e^ > ef) or a hole (e.; < Ef), 
respectively, £f being the Fermi energy. 

To lowest order in a series-expansion of the self-energy 
in terms of the frequency-dependent screened interaction 
W(r, r'; w), the self-energy is obtained by integrating the 
product of the interacting Green function G(r, r'; Si — lo) 
and the screened interaction W(r, r'; u>), and is therefore 
called the GW self-energy. If one further replaces the 
interacting Green function by its noninteracting coun- 
terpart G°(r,r';£i — uj), one finds the G°W self-energy 
and from Eq. {T]) the following expression for the G°W 
reciprocal lifetime: 

rr 1 = T 2]T /dr f dr>*(r)^(r') 
/ 
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x ImW(r,r';\e t -e f \)Mr')<t>f(r), (2) 

where the sum is extended over a complete set of single- 
particle states (j>f(r) of energy e/ (sf < £/ < £j or < 
£/ < £f)- Equation © exactly coincides with the result 
one would obtain from the lowest-order probability per 
unit time for an excited electron or hole in an initial state 
<j>i(r) of energy Si to be scattered into the state </>/(r) 
of energy Ef by exciting a Fermi system of interacting 
electrons from its many-particle ground state to some 
many-particle excited stated 

The interaction W-^r, r'; ui) entering Eq. ((2|) can be rig- 
urously expressed as follows 



W{r,r';uj) 



v(r, r') + J dri / dr 2 u(r,ri 



x x(ri,r 2 ;o;)u(r 2 ,r / ) ) 



(3) 



v(r, r') representing the bare Coulomb interaction and 
x(r,r';tj) being the time-ordered density-response func- 
tion of the many-electron system, which for the pos- 
itive frequencies (u> > 0) entering Eq. {2} coincides 
with the retarded density-response function of linear- 
response theory. In the framework of time-dependent 
density- functional theory (TDDFT), 35 the exact retarded 
density-response function is obtained by solving the fol- 
lowing integral equation)^ 



X(r, r';u>) 



X°(r,r';w)+ J dr a / dr 2 X°(r,ri; 



x {»(n,r 2 ) + r c [n ](ri,r 2 ;w)}x(r 2 ,r';cj), (4) 

where % (r, t';uj) denotes the density-response function 
of noninteracting Kohn-Sham electrons, i.e., indepen- 
dent electrons moving in the effective Kohn-Sham poten- 
tial of density- functional theory (DFT). The frequency- 
dependent XC kernel / xc [no](r, t'uj) is the functional 
derivative of the frequency-dependent XC potential 
V xc {n](r,uj) of TDDFT, to be evaluated at n (r). In 
the RPA, f xc [n ](r, r'; to) is set equal to zero and Eq. @ 
yields the so-called G°W° (or G°W-KPA) reciprocal life- 
time. 

The xc kernel f xc [no](r, r'; uj), which is absent in the 
RPA, accounts for the presence of an XC hole associ- 
ated to all screening electrons in the Fermi sea. Hence, 
one might be tempted to conclude that the full G°W 
approximation [with the formally exact screened inter- 
action W of Eq. |3])] should be a better approximation 
than its G°W° counterpart [with the screened interac- 
tion W evaluated in the RPA]. However, the XC hole 
associated to the excited hot electron is still absent in 
the G°W approximation. Therefore, if one goes beyond 
RPA in the description of W, one should also go beyond 
the G°W approximation in the expansion of the electron 
self-energy in powers of W. By including XC effects both 
beyond RPA in the description of W and beyond G°W in 
the description of the self-energ y 26 ' 27 the so-called GWT 
approximation yields a lifetime broadening that is of the 



G°W form [see Eq. ©], but with the actual screened in- 
teraction W(r, r'; uj) of Eq. replaced by a new effective 
screened interaction 



W{r,r';uj) =v(r,r')+ / dr x / dr 2 {u(r,ri) 

+ r c [n ](r,ri;w)} x(ri, r 2 ; uj) v{r 2 , r'), (5) 

which includes all powers in W beyond the G°W approx- 
imation. 

A. Bounded electron gas 

In the case of a bounded electron gas that is transla- 
tionally invariant in two directions, such as the jellium 
surface or the physically motivated model surface de- 
scribed above, the single-particle states entering Eq. @ 
are of the form 



<£k,i(r) = (j>i{z)e 



ik-ru 



with energies 



£k,i 



k 2 /2m. l 



(6) 



(7) 



k being a wave vector parallel to the surface and m, 
denoting the effective mass in the plane of the surfaced 
Introducing Eqs. ([6]) and (JT]) into Eq. (fT]), one finds 
the following expression for the reciprocal lifetime of a 
quasiparticle (electron or hole) that has been added in 
the single-particle state 0k,i(r) of energy e^: 



'k,t 



j = T 2 dz dzV*(z)ImE(z,z';k, £k!l )^(^), (8) 



where z'\ k, £k,i) represents the two-dimensional 
(2D) Fourier transform of the self-energy S(r, r'; £k,i)- 



1. G°W approximation 

Using the single-particle wave functions and energies of 
Eqs. © and ©, the G°W reciprocal lifetime of Eq. © 
yields 

T k,l = T 2 E/ (02 / d2 / dz>rdz)<f>}(z') 



x lmW(z,z';q,uj)<j) i (z')(tj f (z), 



(9) 



where uj = [(e, + k 2 /2m,i) — (sf + q 2 /2rrif)\, k and 
q represent wave vectors parallel to the surface, and 
W(z, z'\ k, ui) denotes the 2D Fourier transform of the 
screened interaction W(r,r';uj) of Eq. Q, i.e., 



W(z, z'\ k, uj) — v(z, z';k) + / dzi / dz 2 v(z, z\\ k) 
x x(zi,z 2 ;k,uj)v(z2,z';k), (10) 
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v(z, z';k) and %(z,z';k,w) being 2D Fourier transforms 
of the bare Coulomb interaction and the density-response 
function of Eq. ((U, respectively. 

In the G°W° (or G°W-RPA) approximation, the re- 
ciprocal lifetime is also given by Eqs. ([9]) and (jlOp . but 
with the XC kernel f xc [no] (r, r'; uj) entering Eq. (j4|) set 
equal to zero. 

2. GWT approximation 

Using the single-particle wave functions and energies 
of Eqs. (J6]) and ((7|), the GWT reciprocal lifetime is also 
found to be given by Eq. ©, but with W(z,z';k., uj) 
replaced by the 2D Fourier transform of the effective 
screened interaction W(r,r';uj) of Eq. ©, i.e.: 

x ImW(z,z';c l ,Lj)4> i (z')4> f (z), (11) 

where 

W(z, z'; k, ui) = v(z, z'\ k) + / c\z\ J d.Z2 {v(z, zi;k) 

+ f xc [n ](z,zi;\s.,uj)} x(*i, z 2] k, uj) v(z 2 , z'\ k), (12) 

/ xc [no](2;, zi; k, uj) being the 2D Fourier transform of the 
XC kernel f xc [n ](r, r'; ui). 

In the GW°r approximation, the reciprocal lifetime 
is also given by Eqs. (fTTj) and (fT2|) , thereby with full 
inclusion of the XC kernel entering Eq. (fl~2|) . but with 
the XC kernel entering Eq. {4| set equal to zero. 

Hence, we note that both G°W and GWT recipro- 
cal lifetimes [Eqs. ((9|) and (fTTj) ] can be calculated from 
the knowledge of two basic ingredients: (i) single-particle 
wave functions and energies of the form of Eqs. © 
and |(7J), which arc also basic quantities in the eval- 
uation of the noninteracting density-response function 
X°(z,z';k,uj), and (ii) the XC kernel f xc [n Q ](z, z'\ k, uj). 

B. Single-particle wave functions and energies 

For the description of the noninteracting density- 
response function x°(z,z';k, uj) [and, therefore, the 
screened interaction W(z, z';k,uj) and the effective 
screened interaction W(z,z';k,u) of Eqs. (fT0|) and (fT2|) . 
respectively] single-particle wave functions and energies 
can safely be taken to be the eigenvalues and eigen- 
functions of a jellium self-consistent Kohn-Sham hamil- 
tonian.i£ Nevertheless, the actual band structure of sp 
electrons near the surface of noble metals calls for a 
more realistic description of the single-particle wave func- 
tions [<pi(z) and <fif( z )] an d energies [e* and £/] entering 
Eqs. © and fTTj) . 

Hence, in the calculations presented in this paper all 
the single-particle wave functions and energies (those en- 
tering Eqs. §§§ and (jlip and also those involved in the 



evaluation of the noninteracting density-response func- 
tion) are taken to be the cigenfunctions and eigenvalues 
of a physically motivated ID model hamiltonian that ac- 
curately reproduces the projected band gap and the bind- 
ing energy of the Shockley and the first image stated 



C. The XC kernel f xc [no] (z,z';k,u>) 

In order to investigate the impact of strong variations 
of the electron density induced near the surface, and 
because the excitation energies of interest are typically 
small (particularly in the case of Shockley holes), we con- 
sider the following adiabatic (u> — 0) approximations of 
the XC kernel f xc [n ][z,z';k,Lj))^. 



1. Adiabatic local-density approximation (ALDA) 

If one assumes that dynamic electron-density fluctua- 
tions are slowly varying in all directions, the XC kernel 
f xc [no)( z , z '\ k, ui) is easily found to be given by the fol- 
lowing expression: 38 

f xc [n ](z, z'- k, uj) = f xc (n (z); k 3D = 0, u = 0) S(z - z'). 

(13) 

Here, f xc (no(z);k 3D ,u) is the 3D Fourier transform of 
the XC kernel of a homogeneous electron gas of density 
no(z), which in the limit as k 3D — > and lu — > is known 
to be the second derivative of the XC energy e xc (n) per 
particle of a homogeneous electron gas, to be evaluated 
at the local density no (z). We use the Perdew-Wang 
parametrization 3 ^ of the difussion Monte Carlo (DMC) 
XC energy e xc reported by Ceperley and Alder—. 



2. Refined ALDA 

A more accurate description of short-range XC ef- 
fects can be carried out by replacing the local XC ker- 
nel f xc {n (z);k 3D = 0,cj = 0) entering Eq. (UJ) by a 
more accurate still adiabatic but momentum-dependent 
XC kernel f xc (no(z),k 3D = k,u> — 0) (thus only assum- 
ing that the dynamic density fluctuation is slowly varying 
in the direction perpendicular to the surface), i.e, 

f xc [n ](z, z'; k, w) = f xc (n (z); k 3D = k, lj = 0) S(z - z'). 

(14) 

Here we exploit the accurate DMC calculations re- 
ported by Moroni et ai£^ for the static (uj = 0) k 3D - 
dependent nonlocal XC kernel f xc of a homogeneous elec- 
tron gas. A parametrization of this data satisfying the 
well-known small- and large-wavelength asymptotic be- 
haviour was carried out by Corradini et al. (CDOP) 42 . 
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3. Adiabatic nonlocal approximation (ANLDA) 

Here we still neglect the frequency dependence of the 
XC kernel (adiabatic approximation), but now we make 
no assumption on the variation of the dynamic density 
fluctuation and assume that the unperturbed density vari- 
ation [no (z) — no(z')] is small within the short range of 
f xc [no](z, z'\ k, uj). This allows to write 

f xc [n ](z, z'; k, uj) = f xc ([n (z) + n (z')] /2; z, z'; k, u> = 0), 

(15) 



where f xc (n; z, z'; k,ui) represents the 2D Fourier trans- 
form of the XC kernel f xc (n; k, uj) of a homogeneous elec- 
tron gas of density n. An explicit expression for the 
2D Fourier transform of the CDOP parametrization of 
/ xc (n; k, uj = 0) was reported in Ref. [2!J 



r c (n;z,z';k) = - 



4ire 2 C 
kp 



S(z) 



2ne 2 B 
y/gkjr + k 2 



^-y/gk%+k^\z\ 



U3 



2/3- k F z< 



hp 



-0[k 2 F Z 2 /4f3 2 + k 2 /kj 



where C, B, g, a, and (3 are dimensionless functions of the 



Cu(lll)n=0 



Cu(100) n=l 



electron density (see Ref.|42j 
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III. RESULTS AND DISCUSSION 

On the (111) surface of Cu, the n = Shockley state at 
the center of the surface Brillouin zone (k = 0) lies just 
below the Fermi level, with ei — Ef = — 0.445 eV. Bind- 
ing energies of the n = 1 image state on the (111) and 
(100) surfaces of Cu (measured with respect to the vac- 
uum level) are 0.83 and 0.57 meV, respectively. Effective 
masses of the n = 1 image state on Cu(lll) and Cu(100) 
are close to the free-electron mass (rrii = l), 12 while the 
effective mass of the n — Shockley state on Cu(lll) 
is 0.42j 43 ' 44 The probability density of the n = 1 image 
states on Cu(lll) and Cu(100) have a maximum at 2.3 
and 3.8 A, respectively, outside the crystal edge {z = 0), 
which we choose to be located half a lattice spacing be- 
yond the last atomic layer. The n = Shockley state 
wave function in Cu(lll), however, is maximum at the 
crystal edge. 
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FIG. 2: Imaginary part of the screened interaction 
W(z, z'; k, u) and the effective screened interaction 
W(z, z'; k, u), as a function of z = z' and for fixed val- 
ues of k and u) (k = 0.5 A and u = 0.5 eV), in the vicinity 
of the (100) and (111) surfaces of Cu. ALDA calculations of 
Im [W (z, z'; k, oj)] are represented by thick solid lines. RPA 
and ALDA calculations of Im [W(z, z'\ k, ui)] are represented 
by thin solid and dotted lines, respectively. 



A. Screened interaction 

We have carried out calculations of the imaginary part 
of the screened interaction W(z, z'; k, uj) and the effec- 
tive screened interaction W(z, z'\ k, uj) of thin slabs. In 
order to ensure that our slab calculations are a faith- 
ful representation of the actual screened interaction of a 
semiinfinite system, we have used films up to 50 layers of 
atoms and 80 interlayer-spacing vacuum intervals, as in 
the G°W° (G°W-RPA) calculations reported in Refs.QJ 
and[l7l 

The impact of XC effects on the imaginary part of 
the effective screened interaction in the vicinity of the 
(100) and (111) surfaces of Cu is illustrated in Fig. [21 



where ALDA calculations of Imfl/l^z, z'; k, uj)] (with full 
inclusion of XC effects) are compared to calculations of 
Im[W(z,z';k,uj)} with (ALDA) and without (RPA) XC 
effects. Exchange-correlation effects included in the effec- 
tive screened interaction have two sources, as discussed in 
Section Hi A 21 First, there is the reduction of the screen- 
ing due to the presence of an XC hole associated to all 
electrons in the Fermi sea [see Eq. (UJ)], which is included 
in the calculations represented in Fig. [5] by thick solid 
lines and also in the calculations represented by dotted 
lines. Secondly, there is the reduction of the effective 
screened interaction itself due to the XC hole associated 
to each electron [see Eq. (|12p]. which is only included 
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z=z' (a.u.) 

FIG. 3: Imginary part of the effective screened interaction 
W(z, z; k, uS), as a function of z = z' and for fixed values of 
q and u (q = 0.5 A 1 and lj = 0.5 eV), in the vicinity of 
the (111) surface of Cu. ALDA, refined ALDA, and ANLDA 
calculations are represented by thin solid, dashed, and thick 
solid lines, respectively. 

in the calculations represented in Fig. [5] by thick solid 
lines. These contributions have opposite signs and it is 
the latter which dominates. 

Existing GWT calculations of the lifetime broadening 
of image states on Cu(100) and Cu(lll) were performed 
with the ALDA XC kernel that we have used in the cal- 
culations represented in Fig. O The error introduced by 
the use of this local kernel is small in the interior of the 
solid, as the wave vectors involved are small (k < kp). 
However, Fig. [J] shows that the ALDA leads to spurious 
(negative) results for Im[W'(£, z'\ k, us)] near the surface, 
which is due to the presence of small local values of the 
Fermi wave vector (k l p cal < k) in a region where the 
electron density is small. Hence, we have carried out 
refined ALDA and adiabatic nonlocal (ANLDA) calcu- 
lations of Im[Vy(z, z'\ k, lu)) (both with full inclusion of 
XC effects), which have been plotted in Fig.[3l This fig- 
ure clearly shows that the impact of nonlocality on the 
effective screened interaction is large near the surface, 
bringing spurious ALDA calculations (thin solid lines) to 
a more realistic behaviour near the surface (thick solid 
lines). The refined ALDA scheme partially overcomes 
the failure of the ALDA, but a full description of the non- 
locality of XC effects near the surface might be needed 
for a realistic description of the absorption power of solid 
surfaces. 



B. Self-energy 

FigureHexhibits G°W° (G^-RPA), G°W, and GWT 
calculations of the imaginary part of the n = surface- 
state self-energy E(z, z'; k = 0, £k), versus z, in the vicin- 
ity of the (111) surface of Cu, with use (in the case of the 
G°W and GWT approximations) of the adiabatic nonlo- 
cal XC kernel (ANLDA) described in section uTC3l This 
figure shows that as occurs in the case of the screened 





z (a.u.) 

FIG. 4: G°W° (G°W-RPA), G°W, and GWT calculations 
of the imaginary part of the n = surface-state self-energy 
£(z, z'\ k = 0, £k), versus z, in the vicinity of the (111) surface 
of Cu. The solid circle represents the value of z' in ech case. 
GWT calculations (as obtained with the use of our ANLDA 
XC kernel) are represented by dashed lines. G°W (also using 
our ANLDA XC kernel) and G°W° calculations are repre- 
sented by thin and thick solid lines, respectively. ALDA cal- 
culations, which nearly coincide with ANLDA calculations, 
are not plotted in this figure. 



interaction XC effects partially compensate each other, 
leading to an overall effect of no more than 5%. For com- 
parison, we have also used (in the case of the G°W and 
GWT approximations) the ALDA and refined ALDA ker- 
nels described in section Hi C 31 and we have found that 
although the use of these local or semilocal kernels leads 
to spurious results for the screened interaction, our more 
realistic ANLDA kernel yields self-energies that esentially 
coincide with those obtained in the ALDA. 



C. Reciprocal lifetime 

Now we focus on the evaluation of the decay rate (re- 
ciprocal lifetime) of surface-state electrons (and holes) 
at the n = 1 (and n = 0) surface-state band edge 
(k = 0) of the (111) and (100) surfaces of Cu. Cal- 
culations of the noninteracting density-response function 
X°[z, z'\ k, iS) (and, therefore, the reciprocal lifetime) re- 
quire the introduction of complex frequencies of the form 
Lu + b], rj being a positive infinitesimal. Hence, in order to 
ensure that our numerical calculations yield a converged 
value of the reciprocal lifetime, we have calculated r _1 as 
a function of the parameter rj. Fig. [5] represents the re- 
sults we have obtained for the G°W° reciprocal lifetimes 
of Shockley and image states on the (100) and (111) sur- 
faces of Cu, showing that converged results are obtained 
for a sufficiently small value of rj. 

Converged calculations of the reciprocal lifetimes of 
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FIG. 5: G°W° reciprocal lifetimes of Shockley and image 
states on the (100) and (111) surfaces of Cu, as a function of 
the parameter n that accounts for the imaginary part of the 
complex frequencies entering the evaluation of the noninter- 
acting density- response function x°( z > z '\ k, u>). 



son, G a W reciprocal lifetimes are also shown in this ta- 
ble, with use of both the ALDA and the adiabatic non- 
local kernel ANLDA described in section III C 31 Differ- 
ences between our G°W° reciprocal lifetime of the n = 
Shockley state in Cu(lll) (r _1 = 25meV) and those re- 
ported before^ 5 (r _1 = 19meV) are simply due to the 
fact that in our present model we are not accounting for 
the change of the z-dependent surface-state wave func- 
tions 4>i(z) and <pf(z) along the surface-state dispersion 
curve. 

As in the case of the self-energy, the results shown 
in Table [J show that (i) a realistic adiabatic nonlocal 
description of XC effects yields reciprocal lifetimes of 
Shockley and image states that esentially coincide with 
those obtained in the ALDA, and (ii) the overall effect 
of short-range XC is small and GWT reciprocal lifetimes 
are close to their G°W° counterparts. 



TABLE I: G°W°, G°W, and GWT reciprocal lifetimes, in 
linewidth units (meV), of an excited surface-state electron 
(hole) at the n = 1 (n = 0) surface-state band edge (k = 0) 
of the (111) and (100) surfaces of Cu. In the case of the 
G°W and GWY reciprocal lifetimes, both ALDA and ANLDA 
exchange-correlation kernels have been considered. 



Surface 


n 


XC kernel 


G°W° 


G°W 


GWT 


Cu(100) 


1 




17.5 








1 


ALDA 




24 


17 




1 


ANLDA 




24.5 


17 


Cu(lll) 







25 















30 


24.5 







ANLDA 




30.5 


24.5 


Cu(lll) 


1 




29 








1 


ALDA 




42.8 


28.5 




1 


ANLDA 




43 


28 



Shockley and image states on the (100) and (111) sur- 
faces of Cu are exhibited in Table |TJ This table shows: 
(i) G°W° results, which reproduce previous calcula- 
tion o 16 ' 17 , (ii) ALDA GWT results, which in the case 
of the n = 1 image state on Cu(lll) and Cu(100) repro- 
duce the calculations reported in Ref. Qj] (ALDA GWT 
calculations of the reciprocal lifetime of n = Shockley 
states had not been reported before), and (iii) ANLDA 
GWT calculations, never reported before; for compari- 



IV. SUMMARY AND CONCLUSIONS 

We have carried out extensive calculations of the self- 
energy and lifetime of Shockley and image states on the 
(100) and (111) surfaces of Cu, in the framework of the 
GWT approximation of many-body theory. This approx- 
imation treats on the same footing XC effects between 
pairs of electrons within the Fermi sea (screening elec- 
trons) and between the excited surface-state electron (or 
hole) and the Fermi sea. We have included XC effects 
within TDDFT from the knowledge of an adiabatic non- 
local XC kernel that goes beyond the local-density ap- 
proximation, and we have found that these XC contribu- 
tions (in the screened interaction W and in the expansion 
of the self-energy in terms of W) have opposite signs and 
it is the latter which dominates, leading to GWT recipro- 
cal lifetimes that are only slightly lower than their G°W° 
counterparts. 
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